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CONTRACTOR REPORT 


FREE VIBRATION OF RECTANGULAR PLATES WITH 
A SMALL INITIAL CURVATURE 

SUMMARY 


The method of matched asymptotic expansions is used to solve the transverse 
free vibration of a slightly curved, thin rectangular plate. Analytical results for 
natural frequencies and mode shapes are presented in the limit when the dimensionless 
bending rigidity, e, is small compared with in-plane forces. Results for different 
boundary conditions are obtained when the initial deflection is (1) a polynomial in 
both directions, (2) the product of a polynomial and a trigonometric function, and 
arbitrary . 

For the arbitrary initial deflection case, the Fourier series technique is used to 
define the initial deflection. 

The results obtained show that the natural frequencies of vibration of slightly 
curved plates are coincident with those of perfectly flat, prestressed rectangular 
plates. However, the eigenmodes are very different from those of initially flat pre- 
stressed rectangular plates. The total deflection is found to be the sum of the initial 
deflection, the deflection resulting from the solution of the flat plate problem, and 
the deflection resulting from the static problem. 


I. INTRODUCTION 


In many practical applications of plates, it is almost impossible to start out with 
a perfectly flat rectangular plate. It is usually observed that some initial curvature 
is present in one or two orthogonal directions. The effect of various mathematical 
forms of the initial curvature on the vibration of rectangular plates under different 
boundary conditions is the subject of the present study. This report therefore 
complements the earlier work of Adeniji-Fashola and Oyediran [1] in which the effect 
of a small arbitrary initial curvature on the free vibration of clamped rectangular 
plates was discussed. 

The effect of initial deflections on the deformation of elastic plates is treated 
here using small deflection theory for which superposition is applicable. The forms 
of the initial deflection considered include (1) second order polymonial in two orthog- 
onal directions, (2) a second order polynomial in one direction and a trigonometric 
function in the orthogonal direction, and (3) an arbitrary shape which can be decom- 
posed into Fourier components in two orthogonal directions. Results are presented 
for different boundary conditions. In the literature [2], explicit analytical results 
for the free vibration of clamped, perfectly flat rectangular plates are not known. Even 
for the static problems of curved plates , solutions for only a few special cases exist [ 3] . 


It is easily seen that the normalized equations of linear vibration of curved 
plates are readily amenable to the technique of singular perturbations since a small 
parameter, e, multiplies the highest derivative terms. This method has been used 
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successfully by various researchers [4,5,6] to solve the problem of vibration of pre- 
stressed, flat rectangular plates. These authors neglected the corner regions present 
in the rectangular plates. Indeed, few solutions exist in the literature which include 
the corner problem for a rectangular plate [7]. Neyfeh and coworkers [8,9] have, 
more recently, extended the method to solve the vibration problems of prestressed 
circular plates and cylindrical shells. 

The method of singular perturbations is therefore used in the present study to 
analyze the free vibration problem of a slightly curved rectangular plate under various 
boundary conditions and for various mathematical forms of the initial curvature. The 
effects of corner regions are neglected as are the effects of rotatory and shear 
deformations . 

In the next section, the mathematical problem is formulated and the dimension- 
less form of the governing differential equation is obtained. The method of solution 
is presented in Section 3 while the method is applied to a slightly curved plate with 
various mathematiced forms for the initial curvature in Section 4. The report ends 
in Section 5 with a discussion of the results obtained. 


II. PROBLEM FORMULATION 


We consider the linear, transverse, free vibration of a slightly curved, thin, 
rectangular plate. In the absence of shearing prestress, the equilibrium equation of 
curved plates given by Timoshenko and Woinowski-Krieger [3] in rectangular Cartesian 
coordinates is 


DV^ w! 
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( 2 . 1 ) 


where 


w’(x',y’,T) = WQ'(x',y') + w^*(x’,y’,x) 


( 2 . 2 ) 


The deflection WQ*(x’,y') is termed the initial deflection and Wj^'(x’,y*, x) is the addi- 
tional deflection due to vibration, while w’(x',y*,x) is the total deflection. 

In non-dimensional form, equation (2.1) becomes 
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L is the characteristic length used in rendering w', x', and y* in equation (2.1) 
dimensionless. Time is rendered dimensionless by using the characteristic frequency 
0 ) so that 


t = u)T and N^/yto^ = 1 (2.5) 

and y is the mass per unit area of the plate. 

In this work, we shall examine three different cases of initial curvature, viz: 

(1) When w^ is a second order polynomial in both the x and y directions. 

(2) When w^ is a second order polynomial in one direction and sinusoidal in 
the other direction. 

(3) When w^ is arbitrary. 

In this last case, the arbitrary w^ can then be expressed by a double Fourier series. 
If we make the substitution 


Wj^(x,y,t) = u(x,y,t) + fi(x,y) 


( 2 . 6 ) 


in equation (2.3), we obtain two different equations as follows: 
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( 2 . 8 ) 


Equation (2.7) is the governing differential equation describing the dynamic response 
of an initially “flat, normally- prestressed rectangular plate, while equation (2.8) cor- 
responds to the static problem of a similar plate but having an initial curvature, with 
the initial curvature terms on the right-hand side, rendering the differential equation 
inhomogeneou s . 

To complete the problem formulation, we specify the boundary conditions. For 
definiteness , We require that 
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w(0,y) = (0,y) = 0 


w(l,y) = ^ (l,y) = 0 

w(x.O) = 1^ (x,0) = 0 
ay 


(2.9) 


w(x,b) = ||- (x,b) = 0 


0 < b < 1 


for a clamped plate. It is to be noted that the characteristic length, L, has been 
taken equal to the plate length in the x direction. 

Alternatively, the conditions of equation (2.9) can be written as w = 8w/9n = 0 
on all boundaries while the boundary conditions for a fuUy- hinged rectangular plate 
takes the form 


w 


= 9^w/3n^ = 0 


( 2 . 10 ) 


on all boundaries . Here 9 / 9 n denotes the derivative normal to the boundary . 

Equations (2.7) and (2.8) are solved in conjunction with appropriate boundary 
conditions using the method of singular perturbations. The natural frequencies and 
the flat plate component of total deflection are obtained from equation (2.7) while the 
component of the total deflection resulting from the static conditions are obtained 
from equation (2.8). 


III. METHOD OF SOLUTION 


The asymptotic method of singular perturbations is used here to analyze the 
linear, transverse free vibration of a highly-prestressed, initially -curved plate with 
rectangular edges. For infinitesimal transverse vibrations, the effect of mid-plane 
stretching on the in-plane loads can be neglected. With this technique, explicit 
analytical results are obtained in the limit when the applied loads are very large 
compared with D/L^. 

In the limit as e ^ 0^, a state of membrane prevails over the plate, except in 
a thin layer adjacent to the boundaries where higher derivatives become important. 

Under time-independent boundary conditions, we seek for a separable solution 
of equation (2.7) of the form 


u(x,y,t;e) = u(x,y;e) x(t;e) 


(3.1) 
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describing standing waves. Equation (3.1) inserted into equation (2.7) leads to 


2„4 „ 2 9^u „ 2 a^u ,2 _ 

e V u - 3]^ 2 " ^2 2 ' ^ ^ 0 


(3.2) 


ax 


ay 


as the dynamic problem governing the transverse vibration of the initially- curved 
plate. 

The governing differential equations (g.d.e) for the vibration of curved plates 
can thus be written as 
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subject to appropriate boundary conditions. 


Outer Solution 


form 


We seek an approximate solution, otherwise known as the outer solution, in the 


u® = Uq® + eUj^® + e^U 2 ® + O(e^) 


^2 , 2 . , 2 . 2. 2 . n/- 3s 
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(3.4) 


fi® = JIq® + efij^® + £^^2® + O(e^) 


The g.d.e. can then be written as 
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The subscripts denote the order in e, while superscript 0 denotes the outer solutions. 
Variables with negative indices (subscripts) are discarded. Leading and higher order 
equations are given by the recurrence relations equations (3.5) when v takes the 
values 0, 1, 2, and so on. It is pertinent to point out that the order of the g.d.e. 
has been reduced by two from fourth-order equations (3.3) to second-order equations 
(3.5). 


The solutions of equations (3.5) will not, in general, satisfy the prescribed 
boundary conditions. Thus, the solutions are not valid near the plate boundaries. 
Therefore, a thin layer where plate displacements change very rapidly from a membrane 
type [see equations (3.5)] problem to a bending type problem must exist near the 
boundaries. These solutions, valid in the thin layer adjacent to the boundaries, are 
called inner solutions and are matched with the outer solutions using the Van Dyke 
[11] matching principle. 


Inner Solution 


Near the boundaries, where the fourth-order and second-order derivatives are 
of identical orders of magnitude in e, new coordinates (stretched regular coordinates) 
are defined. For example, near x = 0, the new coordinate takes the form 


X = x/e 


(3.6) 


while we seek an expansion of the form 
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where, as before, subscripts denote the order in e and superscript i denotes inner 
solution . 

The g.d.e. then take the forms 


V 1 V 


= -2 u 


XXXX 


XX 


v-2 

XXyy 


“'v-2 

yy 


- u 


v-4 

yyyy 


(A, 


0=0 


2 i 
u 

V- o- 


and 


_y 

XXXX 


- 6 / = -2 


’ V 

XX 


v-2 
XXyy 


+ 3o '!> 


v-2 

yy 


- 

V- 


4 

yyyy 


-(^1= 


+ 3< 


9(eX)' 


3" 

ay* 


)"o 


(eX.y) 


(3.8) 


Expansions similar to equations (3.7) and (3.8) are also written near x = 1, y = 0, 
and y = b, respectively. 

The procedure for obtaining the inner and outer solutions to various orders is 
as follows: If we consider the static problem of equation (3.3), the leading order 

outer solution, is first obtained. The leading order inner solution, is then 

obtained from equation (3.8) and to satisfy all the boundary conditions. These two 
solutions are subsequently matched. The procedure for matching is simple: we write 

the outer solution in inner variables and expand to 1 term while the inner solution is 
written in outer variables and expanded also to 1 term. These two expansions are 

then equated (matched) to determine the unknowns. Next, we solve for and 
carry out a 2-1 matching. The solution of is then sought, with a 2-2 matching 
providing all the unknowns of Next, the inner solution, is determined 

from a 3- 2 matching and , finally , the problem is completed to 0( e^) by the determina- 
tion of from 3-3 matching. 

In the next section, solutions to equations (3.5) and (3.8) are obtained and 
these are demonstrated for various mathematical forms of the initial curvature. 


IV. APPLICATION TO VARIOUS CASES 


The method of singular perturbations (otherwise known as the method of matched 
asymptotic expansions, MMAE) is used in this section to solve the problem of the linear 
vibration of curved rectangular plates. Essentially, the problem has been reduced to 
that of solving two auxiliary problems, equations (3.3), subject to appropriate 
boundary conditions. As stated earlier, equation (3.3a) is identical to the equation 
governing the linear vibration of a perfectly flat, prestressed rectangular membrane. 

We shall henceforth refer to this as Problem I or the dynamic problem. Equation (3.3b) 
is Problem II or the static problem. The initial curvature effect is included in Prob- 
lem II . 

In this section, various mathematical forms of the initial curvature will also be 
examined. 


4.1 The Dynamic Problem 


The governing differential equation for a full- clamped rectangular plate under- 
going sinusoidal vibration is given by 


2o4 

e V u - 


3x 3y 


subject to 
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u(0,y) = 1^ (0,y) = 0 


“ fx ® ’ 

0 < b ^ 1 


(4.1) 


u(x,b) = -|y (x,b) = 0 . 


Schneider {4] was the first to provide an asymptotic solution to the problem 

2 2 

defined by equation (4.1) for 3^ = 3o = 1- The solution he provided is good to 

^ ^ 2 2 2 
0(e). This solution was later generalized for variable 3j^ and and to 0( e ) by 

Hutten and Olunloyo [5]. More recently, Oyediran and Gbadeyan [6] used the same 
asymptotic method to investigate the vibration of a prestressed, rectangular thin 
plate exhibiting natural material orthotropy. Various boundary conditions were con- 
sidered . 

It is pertinent to point out that a similar problem has been solved by Nayfeh 
et al. [8] for circular as well as near-circular and annular prestressed plates. 

We shall merely infer the result to Problem I from the literature. The uniformly 
valid eigenfunction is given [4,5,6] as 


u(x,y) = Aq sin (n-nx) sin (miry/b) 

+ e Aq {(nir/3j^)(2x-l) cos (mix) sin (m7iy/b) 

2 

+ (mir/32l> ) (2y-b) sin (mix) cos (miiy/b)} 

+ Aq {(2n^TT^/3j^^) x (1-x) cos (mix) cos (miry/b) 

+ (nir/3j^^32) ( 232 b- 3 j^) (2x-l) cos (mix) sin (mTiy/b) 

+ (mTi/3j^32^b^) (23j^-32b) (2y-b) sin (mrx) cos (mTiy/b) 

+ (2mV/3o^^) y (b-y) sin (mix) sin (mTiy/b)} 

^ (4.2) 

(continued) 
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+ e {(bj^/6j^) exp(-gj^ x/e) + (bj^/3^) exp(- 3j^(l-x)/e) 

+ exp(-32 y/e) + (- ^ 

+ e {(b2/3j^) exp(-3j^ x/e) + (b2/3j^) exp (- 3j^( l-x)/e) 

+ (^ 2 /^ 2 ^ exp(-32 y/e) + exp(-32 (b-y)/e)} 


+ O(e^) 


(4.2) 


while 


A^(e) = 3^^ TT^ + ^2 + e {4 3^ TT^ + 4 32 ir^/b^} 


+ e^ {(Tr'*/b^) (n^b^ + m^)^ + (12 ir^/b^) (n^ b'* + m^)} + O(e^) (4. 


3) 


where 


bj^ = Aq n IT sin (imry/b) 


n +1 

bj = (-1)"^^^ b^ 


fl = (Ap imr/b) sin (mrx) 
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> (4.4) 
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We shall now consider the static problem (Problem II) for various mathematical 
forms of the initial curvature. 


4,2 The Static Problem 

For a fully clamped plate, the governing equations for the static problem take 
the form 


o 4 

e V - 3 


2 3^ 


3x 
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3y 
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2 ^ '^0 
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1 ^2 
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(4.6) 


Subject to 


«(0,y) = II (0,y) = 0 ; 


fi(l,y) = If (l,y) = 0 


0 < b < 1 


(4.7) 


Qix,0) = If (x.O) = 0 


fUx,b) = If (x,b) = 0 

Here, we shall consider various mathematical forms for w^(x,y) and systematically 

2 ^ 
obtain f2(x,y;e) to 0( e ). 

4.3 Initial Deflection Form — Case 1 

— When the initial curvature is second-order polynomial in both x and y 
directions 

For this case, we let the initial curvature take the form 

Wg(x,y) = d^j^xy (1-x) (b-y) (4.8) 

where is the initial amplitude. 
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4.3-1 Fully Clamped Plate [C-C-C-C] 

It can be seen from equation (4.8) that the initial curvature, w^, vanishes on 

all the boundaries as it must. However, there are discontinuities in the slopes of 
along the boundaries; that is equation (4.8) does not satisfy the relation ^ 


9Wq 9Wq aw. aw„ 

"a^ (0,y) = (I.y) = (x,0) = ^ (x,b) = o . 
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(4.9) 


form 


It is easy to see that equations (4.6) and (4.8) can be combined to take the 


- 32^ ^ = -2 djj 3j^ y(b-y) - 2 b^^ x (1-x) (4.10) 


As outlined previously in Section 3, we write 
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(4.11) 


If we now use equation (4.11) in equation (4.10), we obtain 

2dii[3i^ y(b-y) + x(l-x)] for v = 0 

0 for V = 1 

n ^_2 for V > 2 

(4.12) 

Solving equations (4.12) for provides the outer solution to O(e^). 
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Leading Order Outer Problem ; 

If we set V = 0 in equations (4.12), the resulting differential equation is 
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It can easily be seen that the solution of equation (4.13) takes the form 

A 2 2 2 2 3 

^0 “ ^0 ^2^ + ag xy + X y + 02 x y + Og xy + x y + . . . 

(4.14) 

Clearly, the solution equation (4.14) does not satisfy the boundary condi- 
tions of equations (4.7) and so is not uniformly valid. 

It becomes apparent, after some manipulation, that 

^0 ~ ~ ^3 ~ “4 “ **’ ” ® (4.15) 

while compatibility requires that 

“1 ~ ^11 ’ “2 “ "^^11 ’ “3 ” '^11 (4.16) 

and a^ is yet to be determined. We anticipate that the remaining unknown, a^, will 

be determined when the problem is solved near the boundaries, that is, when the 
inner problem is considered. 


The Inner Problem 

Near the boundaries, we need new coordinates to define the rapidly varying 
deflections accurately. We thus set 


X = x/e , 
X = (l-x)/e 

Y = y/e , 

Y = (b-y)/e 


near x = 0 
near x = 1 
near y = 0 
near y = b 


(4.17) 


Near x = 0 and x = 1, the differential equation for the inner problem takes 
the form 
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Near y = 0 and y = b, the inner problem is obtained by using equation (4.17) 
and following the same procedure, as 
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(4.19) 


Leading Order Inner Problem 

If we set V = 0 in equations (4.18) and (4.19), the resulting differential equa- 
tions are 


’^0 ■ ^1^ '^0 " ® X = 0,1 ; (4.20) 

XXXX XX 

and 

i|»Q - 82 ^ ’^'o “ ® y = 0,b . (4.21) 

YYYY YY 


The solutions to equations (4.20) and (4.21) that satisfy the boundary condition 
equations (4.7) can be written as 



Ki(y) 

K2(y) 

Kg(x) 


' K^(X) 


[X + (l/3i) 
[X + (l/3i) 
[Y + ( 1 / 32 ) 
[Y + ( 1 / 32 ) 


(exp(- 3iX)-l)] , 

(exp(-3iX)-l)] 
(exp(-32Y)-l)l , 

(exp(-32Y)-l)] , 


near x 
near x 
near y 
near y 


= 0 , 

= 1 . 

= 0 , 

= b . 


(4.22) 
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where functions K^(y), K 2 (y), K^Cx), and K^(x) are to be determined from matching. 

Exponentially growing terms have been discarded as unmatchable. Equations (4.14) 
and (4.22), being the solutions to the leading order outer and inner problems, respec- 
tively, are now to be matched. 

From 1-1 matching, it becomes immediately obvious that these two solutions are 
unmatchable. We thus set 


K^(y) = K 2 (y) = K 3 (x) = K^(x) = 0 


(4.23) 


and 


= 0 


(4.24) 


Equation (4.24) is intuitively obvious since a zero value for the leading order inner 
solution is in agreement with the boundary condition of zero deflection on all the 
boundaries . 

To the leading order, the unknown constant, a^, in equation (4.14) remains as 

yet undetermined. It may be determined, however, by considering the first order 
inner problem and, consequently, a 2-1 matching. 


First Order Inner Problem 

The differential equations for are similar to those of iJ;q. Hence, it is easy 
to see that 


Kg(y) [X + 
Kg(y) [X + 
K^(x) [Y + 
Kg(x) [Y + 

The unknowns Kg, 
illustrated below. 


(1/3P (exp(-3^X)-l)] 
(l/3p (exp(-3iX)-l)] 
( 1 / 32 ) (exp(-32Y)-l)] 
(1/32) (exp(-32Y)-l)] 

... Kg and a^ will now 



, near x 
, near x 
, near y 
, near y 

be determined 


= 0 
= 1 
= 0 
= b 


(4.25) 


from 2-1 matching as 


2-1 Matching 

Near x = 0, the 2-term inner solution is 
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(4.26) 


i(X,y) = e 4 ,^ (X,y) + 0 ( e^) 

= KgCy) [eX + (e/ 6 ^) (exp(- 3^X)-1)] + 0 ( e^) 


which can be rewritten in outer variables , x = eX , and expanded to 1 term for x fixed 
and e ->■ O'*" to give 


4'^(x,y) = Kg(y) X + 0 (e) . (4.27) 


Similarly, 1-term outer solution, rewritten in inner variables X = x/e, and expanded 
to 2 terms for X fixed and e O’*" yields 


il°(X,y) = e [agy + ag y^] X + 0 (e“) 
or, rewritten in outer variables, 

^ 2 ^(x,y) = (a^y + y^) x + 0 (e^) (4.28) 


Van Dyke's [11] matching principle states that these two expansions, equations (4.27) 
and (4.28), are equal. Hence, we have that 


K 


gCy) = agy + ag y 


(4.29) 


Similarly, near x = 1, we set 1-x = eX. Using Van Dyke's principle, it is easily 
verified that 


aj = djjb 


and 


(4.30) 


Kg(y) = -d^j^y (b-y) , Kg(y) = Kg(y) 


(4.31) 


The same procedure can be carried out near y = 0 and y = b. It can be readily 
shown that K,^(x) and Kg(x) are determined as 
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= Kg = x(l-x) 


(4.32) 


To complete the first order problem, the outer solution is also required to 0(e). 
First Order Outer Problem 

From equations (4.12), the first order problem takes the form 


3^^ + 32^ oj =0 . 


XX 


yy 


The resulting solution, , is to be matched with 

- Kg(y)/3j^ , near x = 0 

- Kg(y)/3j^ , near x = 1 

- K,^(x)/ 32 ’ y = 0 


- Kg(x)/32 ’ y = b 


(4.33) 


(4.34) 


The requirement that match with equations (4.34) is equivalent to a 2-2 matching. 

Alternatively, the problem defined by equations (4.33) and (4.34) can now be 
rewritten as 


+ 32 ^ =0 . 

XX yy 

Subject to 

= (dj^j^/3j^) y(b-y) ; near x = 0 and x = 1 

x(l-x) ; near y = 0 and y = b 


(4.35a) 


(4.35b) 


We seek a solution to equations (4.35) of the form 
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(4.36) 


~ (<^ 11 /^ 1 ) y(t)-y) + (d^j^b/gg) X(l-x) + $(x,y) 


such that $(x,y) is determined by inserting defied by equation (4.36) into equa- 
tion (4.35a). Thus, the equation for $(x,y) takes the form 


61 ^ $ + <I> + D = 0 

XX ^2 yy 

where 

D = -(2d^^/3^32> + ^2^^ 


(4.37) 


(4.38) 


The form of the solution suggested in equation (4.36) ensures that Q®(x,y) is 
matched along the boundaries provided $(x,y) vanishes on the boundaries. 

It is interesting to note from the symmetric end conditions of equations (4.35b), 

that 


Tx 


•x=l/2 




= 0 


'y=b /2 


(4.39) 


It can easily be shown that 


t(x,y) = S [*!„ sinh(-g^y) . cosh(-j^ y] - 


x) * A,„ oosh(!!^^ x) - (=^) 


(4.40) 


where the constants B^, A 2 jjj, and are determined from a 2-2 matching and 

are given below as 


®n = 


n 

„2 2 ^ 2 

n TT 3 j^ 





(4.41) 

(Continued) 
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m 


*3m b^)] 


'3m „2 2g 2 

m ir ^2 


'2ni 


( mirg„\ 


D = — {(-1)*'-!} and E = ^ {(-1)'"-!} 
n mr m tnir 


(4.41) 

(Concluded) 


This completes the problem to 0(e). 


Second Order Correction 

2 

Since we seek a solution correct to 0(e ), we have to solve the second-order 
problem. This is accomplished by first solving the inner problem and then the outer 
problem . 


Inner Problem 

Near x = 0, the solution to the inner problem of equation (4.18), to second 
order, is 

^ = KgX + d^^y(b-y) x^ + £ (KgX - Kg/ej) - Kg/g^ (4.42) 

where Kg(y) is given by equation (4.31) and Kg(y) is to be determined from a 2-3 
matching, given below. 

2- 3 Matching 

The 2-term outer solution 


fi®(x,y) = fig^(x,y) + e fij^®(x,y) + O(e^) , 


rewritten in inner variables near x = 0, say, and expanded to 3 terms for X fixed 
and e -> O'*" is to be matched with the appropriate terms of the inner solution given 
by equation (4.42). It is easily verified that, equating the ex terms, yields 


Kg(y) = dj^j^b/32 ^ “l^ ®n “l^ “ ^2m “2 (miry/b) 

(4.43) 
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where 


= (mr3j^/32) . “2 “ (mTr32/b3j) • (4.44) 

The same procedure can be carried out near x = 1, y = 0, and y = b, respec- 
tively. It is easy to see that the second order outer solution, ^ 2 ^* have to 
match with -Kg/3j^ near x = 0. This is equivalent to a 3-3 matching. 

Outer Problem 

For a 3-3 matching, it is obvious that ^2^ must be obtained from the differential 
equation 


“2“ '“2“ = ‘“0“ 


(4.45) 


XX 


yy 


and must then be matched with the inner solution as 


*2^(0, y) = t2^0,y) = -(d^j^b/33^32) ‘ (nir/3p sinh a^^y + (cosh a^y-1)] 


(A2a2/32) sin (mtry/b) 


iJ2°(l»y) = t2^1.y) = -(d^j^b/3j^p2^ (nir/e) [Aj^^ sinh a^^y 

+ Bln aj^y-1)] + (I/61) [A2^«2 cosh 02 

+ A2^a2 sinh 02] sin (miry/b) 


a2®(x,0) = ij>2^(x,0) = -(diib/3i32) " (mir/b32) ^^2m “2^ 

+ Ajj^ (cosh o^gX-l)] - (Aijjai/32) sin (mrx) 


(4.46) 

(Continued) 
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fi2°(x,b) = (|<2V,b) = -(d^j^b/B^g2^ (mTr/be2) [A2 j„ sinh a2X 

+ (cosh a 2 X-l)] + (I/B 2 ) cosh cTj^b 

+ oT, sinh cT,b] sin (mrx) . (4.46) 

^ (Concluded) 

To solve this boundary value problem, it is convenient to introduce the decom- 
position 


0 = 6^1) + 6 (2) ^ . (3) 


^ 


(4.47) 


where is a solution of equation (4-45) but does not satisfy equation (4.46). 

We insist here that i >2 ^ satisfy instead the following boundary condition extracted 

from equations (4.46) 




(|.2^^\l.y) = -(d^^b/ 6 j 62 > ®ln^^P 


(|)2^^\x,0) = -(dj^j^b/6j32) 


<P2^^\x,b) = -(dj^b/3j32) ^3m^^^2^ 


(4.48) 


Similarly, the problem for <|> 2 ^°^^ is described by 


a 2* (a) 

2 ^ '*’2 


ax 


2 ^'^2'“' 

2 ^2 2 " ® 

ay 


a = 2,3 


(4.49) 


while the boundary values of (f>^^ ^ (a = 2,3) are derived by insisting that <j> 9 ^ ^ + 

( 2) ( 3) ^ 0 ^ . 

<Pn ‘f*2 ^ must satisfy all the boundary conditions on . Here, we require 

( 2) ^ 
that <t >2 must satisfy the following boundary conditions: 
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(4.50) 


<j>2^^^(0»y) = (rniry/b) 

( 2) 

(j>2 = (02/3^) [A2^ cosh 02 + Agjjj sinh 02] sin (rrnry/b) 

<l'2^^^(x,0) = "(Aj^j^aj^/32) (nirx) 

<l'2^^^(x,b) = (aj^/32) “1^ '*’ ®ln (nirx) 

With little effort, the solution to <^ 2 ^^ obtained as 

4 > 2 ^^^ = ({>2^^' (x,m) sin (rmry/b) + <i>2^^'(ri»y) sin (nirx) 

where 


(j>2 '(x,m) = Fj^ sinh 02X + F2 cosh a2X - K^/a2 


(j>2^^\n,y) 


= Fj sinh ctj^y + F^ cosh a^y ~ Gj^ 


while 


2m^ ir^ 3„ A- . 

j^i = _2_3m ^ | 


Odfib niiT32 ‘^11 


tnir 


3i 


~ ba« ( i K. ^ 

= 2m?l7 > ^^2m “2 “2^ ‘ ^2 

1 ( 1 02 


^2 " ITZ 


^ fi-r-n™! ^2m“2 
5m¥ > ~T[~ 


^ Gj + Gj - cosh oTj^b 
F = 

sinh cTj^b 


(4.51) 


(4.52) 


cosh 02 


(4.53) 

(Continued) 
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^4 = * °2 


2 t ((-1) 


n 

^ y niT 


nir 


«2 2 . 2 „ 2 
n IT ^2 


°2 = 55^ (l-(-l)"] 


°3 = la^ lAln “l*> ■" ®ln “l'’ 


The solution to equation (4.49) for a = 2 can be easily shown to take the form 


= Hj^(y) sin (mrx) + H 2 <x) sin (mry/b) 


It clearly follows from equations (4.54) and (4.49) that 


(|>2 = [Lj^ cosh a^y + L 2 sinh Oj^y] sin (mrx) 


+ [Lg cosh a 2 X + sinh OgX] sin (rnwy/b) 


''In 

h 


H = - 


^2m“2 


2 «2 
L, + — 

^2 “ 1 ^ ^2 


~ _ 2 «2 A2p, a2 

LtM — + 

tanh 02 
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definltionir‘’“’'®“°'’ equation (4.49) for a = 3, we first make the following 



(nirx) cosh a2X dx 


(mrx) sinh S2X dx 


(mTry/b) cosh cTj^y dy 


(imry/b) sinh a^y dy 


(4.57) 


Also, 

^ _ niTT . 

“ ■ b3^ ®2 ^ 3 m 

^2 - ■ b 32 ®2 ^1^ 

(4.58) 

^3 17 ‘Afn «4 ^ Bi„ 63) 

7 - (-1)^ Htt 

4 ®4 ■*■ ®ln 

Then takes the form 

<l>2^^^x,y) = i‘2^^^(n,y) sin (mrx) + ^2^^^x,m) sin (mrry/b) (4.59) 
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where 


<}» 2 ^^^(n.y) = Lg sinh + Lg cosh cTj^y + L,^y sinh a^^y + Lgy cosh a^^y (4.60a) 


7 (3) 


'’(x.m) = Eg sinh a 2 X + Eg cosh S 2 X + E^x sinh a 2 X + Eg cosh a 2 X (4.60b) 


respectively . 

From a 3-3 matching, it can easily be established that 


■[ 


sinh Oj^b 


- Lg coth otj^b - L^b - Lgb coth 


“1^] 


^6 ’ ^7 “ 


n2 .2 g 


1 B L - 


“1 ^2 


2 “In ’ “8 B 


In 


E, 


=[ 


sinh a. 


- (Eg + Eg) coth - Ey 


(4.61) 


®6 ■ ^3 ’ ^7 ■ 


2 . 
m IT 3. 




^7 ^2m 


“2 


2 3m ’ “8 A 


3m 


2 

This completes the static problem to 0(e ) for a clamped, rectangular plate 
whose initial curvature is second-order polynomial in x and y. 

The result for a fully hinged, rectangular plate is now presented and is found 
to check with the exact solution presented by Timoshenko and Woinowski-Krieger [3]. 


4.3-2 Fully Hinged Plate [S-S-S-S] 

For this case, the dynamic problem has been solved [5] and it is found that 
the eigenvalues can be written down exactly as 


^2 _ _2 2 . 2 . _2 2 . 2 . .„4 4 , 4 . „ 2„2 4 . 2 , „4 4 . ,.4 
A -nirgj^ + tnir32 "*■ £(^^ i^b + nmirb + mTr)/b 


(4.62) 


while the dynamic eigensolution takes the form 
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and the static deflection is given by 


o-oOj. 2_0 

n - + e ^2 


(4.64) 


Here, 


^^0 = "^11 Cl-x) (b-y) 


(4.65) 


while ^2^ determined from the equation 


2 


+ 6. 


2 ^ ^2 _ „4 . 0 


3y 


= 7’ n " = -8d 


11 


(4.66) 


subject to £2 2 = 0 on all the boundaries. Therefore, we have that 


^(x,y) - 2 cosh (a^y) + F2 sinh (cTj^y) - L2/ctj^^] sin (mrx) 


+ 2 


L [Fi 

m 


cosh (a2x) + F2 sinh (CI2X) - L2/a2^] sin (irnry/b) 


where 


8 dii 5 

^2 = ■ ^ * C = [1 + (-l)'^ 


^2 =- 


8 djl b ? 
m IT 


, ? = [1 + 


2 


(4.67) 


> (4.68) 


- ^2/0^ . F^ = [Fj^(l - cosh (oj^b))] /sinh (oj^b) , 

’ ^2 ~ ^^1^^ ~ cosh (02))] /sinh (02) . 
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It can be readily verified that 




0 

V 


= 0 


for V = 3,4, .. . 


The static problem for a fully hinged , rectangular plate whose initial curvature 
is second-order polynomial in x and y is thus completed to all orders. 


4.3-3 Partly Clamped and Partly Hinged Plate [C-C-S-S] 

For this case, we shall infer the result to the dynamic problem from the litera- 
ture. This is already given as equation (4.2) above and will not be repeated here. 


The Static Problem 


For a plate that is partly clamped and partly simply- supported on adjacent 
sides, the outer problem to solve takes the form 

2 y(b-y) + x(l-x)] , v = 0 

0 , V = 1 (4.69) 

, V ^ 2 




3x 


+ 6 


2 ^ 


in place of equation (4.12) and subject to 
'l'^(0,y) = t|<^(x,0) = 0 , 

d jp d ]p 

(0,y) = (x,0) = 0 , 

'/'^(i.y) = 4»^(x,b) = 0 

— ? (i,y) = — r (x,b) = 0 
ax*^ ay^ 

Here 

n (x,y) H tp (x,y) 
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> (4.70) 


! 


Following the procedure outlined previously in Section 4.3, the solutions to the 
outer and inner problems are summarized below: 

'I'qCx.y) = -dj^j^ xy(l-x)(b-y) , (4.71) 

'l^l(x,y) = (d^^/3^) y(l-x)(b-y) + (dj^j^/32) x(l-x)(b-y) + (f>^(x,y) (4.72) 

where 4>j^(x,y) satisfies the Poisson equation subject to homogeneous boundary condi- 
tions. With relative ease, it can be shown that 

<^j^(x,y) = cosh (oTj^y) + sinh (a^y) - (G j^/a]^^)(b-y) 

+ (b^ G 2 /n^iT^)} sin (nirx) + {Ej^ cosh (o^gX) + I 2 sinh (ot 2 X) 

- (Gj^ + § 2 ) (1-x)} sin (mTry/b)/b , (4.73) 

where 

Gf = (2dj^j^3j^^/mrg2^) 5]^ . = dj^j^b^/Trg2 

G 2 = 2dj^j^/mr3j^ , G 2 = (2d^^j^b 32^/01^ 

E^ = - (b^G2/mV) . 

(4.74) 

Eg = E^ coth (aj^b) + b^Gg/mV sinh (a^^b) , 

ll = (^1 + » 

Eg = 1^ U - cosh (cT 2 )]/a 2 ^ sinh (oTg) - Gg coth 
Similarly, we have 

4»g(x,y) = tpg^“^(x,y) + 4<g^^\x,y) + il^g^^^x.y) (4.75) 


27 


where 


2 

■ Pq ■ 2^2^ ■ y(Pi+P2y)]/“i^> sin (mrx) 

+ {Fj^ cosh (ct2x) + F2 sinh (a2x) 

- [(^2^ Pq - 2P2) - tt2^ x(P^+P 2^)J/«2^> sin (imry/b) ( 4 . 76 ) 


with 


Pq = (a^2/3^mr) [S^ + 82 b - ( 8 d^^ 6 ^/ 0 ^^ , 


P^ = (a^^/6^mr) (82 - Sgb) , 


P 2 = - (aj2 ^ 11 / 61 ^ ♦ 


Pq = i'^2 ^'^2 ^2 " ^2^^2^ ’ 


P^ = («2^ b/p2 mir) (§2 - S3) , 






(4.77) 


Fj^, F2, Fj, and F2 remain to be determined. 


( 6 ) 


1^2 (x,y) = (Sg/32) [coth (oj^b) sinh (oj^y) - cosh (oj^y)] sin (mrx) 


+ (Sg/ 3 j^) [coth (02) sinh (a2x) - cosh (0^2^)! sin (miry/b) 


(4.78) 
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with 


Sg = 2 [a^ + G^/aj^] 

Sg = 2 [Sj §2 + Gj/Sj^] 


and 


(4.79) 


(y) 


'^2 t(Nj^ + NgY) cosh (oj^y) + (N 2 + N^y) sinh (otj^y)] sin (ntrx) 


+ [(N^ + Ngx) cosh (a 2 x) + (Ng + N^x) sinh (St 2 x)] sin (imry/b) 

(4.80) 


where 


N 3 = Ei/3 , 

Ni = -(S^ 0 j^ + Sg 62)732 

N 2 = - [N^b + (Nj^ + Ngb) coth (oj^b)] 

= 2 m TT i^/b^ , Sg = 2 tn IT 12 /b^ 

= 2 n ir b^ G 2 /m^ , 

®2 = ‘^11^^2 ■ 2 n TT G^/a^2 ^ 

®3 ~ ’ ®4 ~ 2 n IT Ej^ , Sg = 2 n ir E 2 

= 2 m IT 02 /b^ ir^ , 

®2 " - 2 m ir G^/b^ 02 ^ 
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^3 “ ‘ 


and 


Ni = - (S^e^ + S 202 )/ 3 i 


^2 = ■ ^^4 + (Nj + Ng) coth (S2)] 


N3 = 012 S2 b/2 n Ti 32 


N^=a 2 Sj^b/ 2 mTr 34 


( 4 . 81 ) 

(Concluded) 


J 


We now determine the remaining constants, F2 as follows. We set 


and 


Vj = (-<Sj + Sj) Cl + (- 1)"^^ (§2 - S3) 


- Sg (2?j^ + (-1)™ n^ Tr^)/n^tr^}/ 82 


V2 = {-(S^ + 82b) + (-l)'"'^^ (S2 - Sgb) 


- S3 (2?^ + (-1)"” TT^) b^/m^ b/3j^ nnr 


and define 


Fi = V2 + (Si^ P3 - 2P2)/ffi^ 




( 4 . 82 ) 


F2 = -F^ coth (Jjb) + [(a^^O ' ^^ 2 ^ " A 4 . 83 ) 


^1 = ^2 + Po - • 


F2 = -F^ coth (J2) + [«2^ Pq ■ ^^2 ■ ^2^^ ^“2^ 
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Finally, the solution to the inner problem to 0(e ) takes the form 


'l'(x,y) = ^ 


Tj^(y) X + y(b-y) + e[Tj^(y) x - Tj^(y)/3^] - Tj^(y)/gj^ 

near x = 0 

T2(y) (1-x) + y(b-y) (1-x)^ + e[T2(y) (1-x)] , near x = 1 


T,(x) y + d^^ x(l-x) y2 + eCT^fx) y - T,(x)/3<^] - ^ T,,(x)/3, , 


3v^// 1^2 

near y = 0 


T^(x) (b-y) + dj^^ x(l-x) (b-y)'^ + e[T^(x) (b-y)] , near y = b 


where 


(4.84) 


Tj^(y) = T2<y) = -d^^ y(b-y) , 

TgCx) = T^(x) = -d^j^ x(l-x) , 

Tj^(y) = Sj^ + S2(b-y) + Sgy(b-y) + cosh (aj^y) + Sg sinh (oj^y) 

+ Sg sin (mTiy/b) , 

T2(y) = (-1)"'^^ + S^(b-y) + S3 y(b-y) + (-1)*''^^ [S^ cosh (cTj^y) 

+Sg sinh (cTj^y)] + Sg sin (miry/b) , 


T3(x) = Sj^ + S2(l-x) + S3 x(l-x) + S^ cosh (a2X) + Sg sinh (a2x) 

+ Sg sinh (mrx) , 

f^Cx) = S^ + S^(l-x) - S3 x(l-x) + (-1)“'^^ [84 cosh (ogx) 

+ Sg sinh (a2x)] + Sg sin (mrx) 


> (4.85) 


This completes the static problem to O(e^) for a rectangular plate that is partly 
clamped and partly hinged on adjacent sides and whose initial curvature is second- 
order polynomial in x and y directions. 
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4.4 Initial Deflection Form — Case 2 


When the initial curvature is polynomial in one direction and sinusoidal in the 

other : 

Here, we set the initial curvature to take the form 

Wg(x,y) = dj^j x(l-x) sin (mity/b) (4.86) 

with dj^j^ being the amplitude. 


4.4-1 Fully Clamped Plate [C-C-C-C] 

As for the previous case, satisfies the zero deflection condition on all the 

edges while the condition of zero slope is not met. For a fully clamped, curved plate, 
the problem can be decomposed into two, viz: the dynamic problem and the static 
problem. The former is the same as was presented in equations (4.2) through (4.5), 
while the latter problem takes the same form as equations (4.6) and (4.7) with equa- 
tion (4.8) being replaced by equation (4.86). 

It is easily verified that the leading order solution takes the form 


Q 



11 


x(l-x) sin (miry/b) 


(4.87) 


while the inner solution, to 0(e), can be written as 


V = 


e Ujj(y) (X + ( 1 / 6 ^) (exp(-3jX) - 1 )) 

e Uj^(y) (X + (l/3j) (exp(-gj^X) - 1 )) 

e U2(x) (Y + (1/62) (exp(-32Y) - 1 )) 

e Ug(x) (Y + (1/62) (exp(-32Y) - 1 )) 


, near x = 0 


, near x = 1 


, near y = 0 


, near y = b , 


(4.88) 


where 


UQ(y) = Uj(y) = - dj^j sin (miry/b) 


U2(x) = Ug(x) = - (dj^j^mir/b) x(l-x) 


(4.88a) 
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The outer problem to 0(e) is thus identical to that of Case 1, equation (4.35a), but 
with the following boundary conditions: 


S^j®(0,y) = (dj^j^/pj^) sin (miiy/b) 
aj^®(l,y) = (d^j^/g^) sin (miry/b) 
flj^®(x,0) = X (1-x) 

J2j^®(x,b) = (dj^j^/ggb) X (1-x) 



The solution to equations (4.35a) and (4.89) is easily obtained as 


fij^^(x,y) = (dj^j^/gj^) sin (miry/b) + (dj^j^ mii/g 2 b) x (1-x) + l'^(x,y 


such that $ is determined from 


9 9 

^1 n + h ^ 

XX yy 


. -( 


2dj^jmTTgj^^ g2^djj^m\^ 

• + R sin ( 


^ 2 ^ 




T>) = 


subject to # = 0 on all the boundaries. 

The solution to equation (4.91) is obtained as 

$ cosh (a 2 x) + Kj ^2 ^2^ (miry/b) 

+ [Kj^j cosh (o^y) + sin (oTj^y) + Ig] sin (mrx) 


+ sin (miry/b) sin (mrx) 


where 


^2 = 


2 dj^ [l-(-l)“] b^ g^2 


m3 .2 « 3 


(4.89) 


(4.90) 


(4.91) 


(4.92) 


(4.93) 

(Continued) 
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2 ^11 ^2 
b 3g IT 


- (-l)"]/nir 


(4.93) 

(Concluded) 


K = 

V 


62^ m2 


The unknowns ..., are determined from a 2-2 matching and are given 


*^11 ^2 ’ ^12 


*^13 ^3 ’ ^14 


N2 (cosh 02 - 1) 


sinh a„ 


(4.94) 


Ig (cosh (oTj^b) - 1) 
sin (Sj^b) 


The problem is thus completed to 0(e). 


Second Order Correction 


It is easily shown, from equation (4.12) and 3-2 matching, that must 
satisfy the equation ^ 


2 *^“2° 


2 ’^“2" 


! 2 ^ 2 "2 _ „4 . 0 

1 . 2 ^ ^2 ^2 “ ^ ^0 

3x 3y^ " 


(4.95) 


‘=2°«'’y> = ^ ♦ "3 "’) - ^ <*12 ®2 * nu) sin (5!^) 


[Kj^j cosh (otj^y) + sinh (oj^y)] 


«2'’(i.y) = - i ( 


^ '3 




[Kfi a2 sinh a2 + Kj^2 “2 “2 


+ (-1) rnr] sin (5^) + (-1)” [Kj 3 cosh (a^y) + sinh (a^y)] 
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~ ■ $2 \ ^2 ■ ^2 \^14 “l 


niTT , „ miT 


^ sin (niT 


X) 


^ (“aX) + K,o sinh (a,x)] 


‘12 


0 


”(x,b) = 


_ (- 1 ) 


m 




™ It ^^11 ^ ^12 ^“ 2 ^^^ 

2 


+ ^^13 “l (oj^b) + a-^ cosh (cTj^b) + (- 1 )™ K^ 5~] sin (mrx) 


(4.96) 

(Concluded) 


The solution to this problem is analogous to that of the problem defined by 
equations (4.45) and (4.46) and can thus be inferred directly from equations (4.47) 
through (4.61). Thus, this completes the static problem to 0 (e 2 ) for a clamped, 
rectangular plate whose initial curvature is second-order polynomial in one direction 
and sinusoidal in an orthogonal direction. 


4.4-2 FuUy Hinged Plate [S-S-S-S] 


For a fully hinged plate having its initial curvature given by equation (4.86), 
the solution to the associated dynamic problem is the same as given in equations (4.62) 
and (4.63). The result for the static problem is now summarized below: 

2 

We have, for the outer solution given to 0(e ), 


4'o(x.y) = - dj^j x(l-x) sin (miry/b) (4.97) 

ij;j(x,y) = 0 (4.98) 

4 » 2 (x,y) = B sin (miry/b) + D x(l-x) sin (miry/b) (4.99) 


where 


B = - A m^ ir^/b^ 62 ^ ^ = < 4 / 82 ^ + 2 8 ^^ A/ 82 ^) 


(4.100) 


with A = The inner solutions all vanish. 
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4.4-3 Partly Clamped and Partly Hinged Plate [C-C-S-S] 


For a plate that is partly clamped and partly simply- supported on adjacent 
sides, the outer problem to solve is given by the differential equation 

2 '^11 [l ^ 2- ^ *<!-='>] sta . 

^ ° V = 0 

0 , V = 1 

V2 ’ ^ - ^ 

(4.101) 


subject to the boundary conditions of equations (4.70). The method of solution is 
the same as was enumerated before and the solution is summarized below; 


= ■ dj^jx(l-x) sin (miry/b) , (4.102) 


i()j^(x,y) = (dj^j^/pj^) (1-x) sin (miry/b) + (dj^j^ mTr/32b2) x(l-x) (b-y) 

+ [Ej^ cosh (a^y) + E2 sinh (oj^y) - (Gj^/ctj^^) (b-y) 

+ (b^/m^ ir^) G2 sin (miry/b)] sin (niix) 

+ [E^ cosh (H2X) + §2 sinh (S2X) - {(G^/a2^) + (G2(l-x)/a2^)}] 

sin (miry/b) . (4.103) 

2 

The unknown constants here are defined later on in the section. The 0(e) 
term, is decomposed into three parts as 

'l'2(x,y) = <|'2^‘^^(x,y) + ij>2^°\x,y) + t|<2^^\x,y) (4.104) 

where 

cosh (cTj^y) + W2 sinh (oj^y)] sin (nirx) + [W^ cosh (012X) 

+ sinh (o2x) + (d^^^ g^^/32^) {(4 ~ ^ ^2 x(l-x)}] sin (miry/b 

(4.105) 


9 9^ 4' 

2 V 

^ 3x^ 


+ e 


O 9^ Ij) 

2 

2 ~~Y 
9y 
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with 


Wi = - (V^ 2 > 

W2 = - coth (cTj^b) 

W3 = -[(S^/ep + (2 (2 ^2 ■ 

= -[W3 cosh + (2 ^2 ~ 3i^)^/sinh (02) . 

\p 2 ^‘^\x,y) = [Nj^ cosh (a^y) + N2 sinh (a^^y) + N3 y cosh (a^^y) 

+ y sinh (a^^y)] sin (mrx) 

+ [Nj^ cosh (ct2x) + N2 sinh (a2x) + N3 x cosh (a2x) 
+ X sinh (012X)] sin (imry/b) 


> 


I 


(4.106) 


1 

i 


with 

N3 = (Hj^ S2/2mrg2^) 

= (otj^ Sj^/2mr3j^) 

Ni = - (S^ 01 + §2 02)/e2 

N 2 = -[N^b + (N^ + N 3 b) coth (oj^b)] 

N 3 = (a 2 §2 b/2mT732) 

= (tt2 Sj^ b/2mnr32) 
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Ni = - (S^ + S 2 




(4.108) 

(Concluded) 


and 


4 ^ 2 ^ ’^^(x.y) = [Fj cosh (Hj^y) + F 2 sinh (cTj^y) - (S 3 /Bj^n 7 r)(b-y)] sin (mrx) 


+ [Fj^ cosh (tt 2 x) + F 2 sinh (ct 2 X) 


(4.109) 


+ {2 P 2 + 0 .^ Pq - 02 * x(P^ - P 2 x)}/a 2 ‘*] sin (imry/b) 


^ 4 . 


with 


P. = (a„^ b/&„ imr) (S. - S.) , 


P^ = (02^ b/82 nnr) (S4 + S3) , 

^2 ~ ^“ 2 ^ ^^^2 ^3 ’ 

F^ = (2 P 2 + ~ ^®3 ^^^^1 ’ 

^2 ~ “2 ^ ^ ^2 ^ “2^ ^2 ^^1 ^2^ ^^^“2^ “2^ 

F^ = {(Sg - S 4 ) q + (-l)"^*^^ (S 4 + S 3 ) + (S 3 /n^A (2 5i 
+ (-1)” ri^ TT^)} / ( 3 „n 7 r) - (S„ b/ 3 -.nTr) 


> (4.110) 


F„ = - F. coth (a.b) 


Also, we have that 


= (2 3 ^^ m d^^/32^ b^n) = (m^ tt d^^/ 3 ^ b^n) , 


= (G^ b/a^^) 


; E 2 = -(G^ h/a^n coth (a^b) , 
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= (2 djj/Bj) 

= (Gj + G 2 )/ 52 ^ 

= 2 n IT Ej^ 


~ ^ *^11 ’ 

^2 “ ^“ 2 ^^ ” ^1 ^“ 2 ^^ 

S 2 = 2 n TT E 2 , 


= n Tr/32b^) - (2 n ir Gj^/a^^) 

= (2nb^ G^Ito^v) + (2 ^ 2 '^) ~ (d^^/3i> + (2 G^ta^) , 

= (d^jin + (-1)” (2 n TT G^/a^^) 

= (dii/3p - (2 0*2^ R2/b) - (2 G2/b a^) , 


> (4.111) 


where 



^2 ” ®1 (“ 2 ^ ^ ®2 (“ 2 ^ 

Finally, the solution to the inner problem complement of equation 
takes the form 

C 2 

bj(y) X + d^^ X sin (miry/b) + e [b 2 (y) x - b^(y)/3j^] 

- b 2 °(y)/ 3 i . 

bj”(y) (1-x) + dj^^(l-x)^ sin (mTry/b) + e b 2 ”(y) (1-x) , 


(4.101), to O(e^), 


near x = 0 ; 

near x = 1 ; 


l'(x,y) 


(x) y + e [f 2 *(x) y - f^*(x)/32] 
fj"(x) (b-y) + e f 2 "(x) (b-y) , 




where 


near y = 0 ; 

near y = b . 
(4.112) 
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bj^*(y) = - sin (tmry/b) 


b^”(y) = b^’(y) , 

fj^*(x) = mit/b) X (1-x) , 

f^”(x) = (-1)“ fj’(x) 


> 


and 

b2*(y) - cosh (ctj^y) + S2 sinh (Sj^y) + (b-y) 

+ sin (miry/b) 

b2”(y) = (-1)”'*’^ [Sj^ cosh (cTj^y) + S2 sinh (ctj^y)] + Sg (b-y) 

+ Sg sin (imry/b) , 

f2*(x) = cosh (a^) + §2 sinh (ot2x) - Sg x(l-x) !> 

+ (1-x) - Sg + Sg sin (niTx) , 

f2*’(x) = (-1)™^^ [Sj^ cosh (ct2x) + S2 sinh (012X)] + Sg x(l-x) 

+ (-l)™"^^ S^(l-x) + (-1)“ Sg - Sg sin (nirx) . 

This completes the static problem to 0(e ). In the next section, we 
the same plate problem with an arbitrary initial deflection. 


(4.113) 


(4.114) 


consider 
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4.5 Initial Deflection Form — Case 3 


- When the initial curvature is of arbitrary form. 

Here, we decompose the initial curvature into its Fourier components, i.e.: 

00 00 

Wjj(x,y) = X) Piun (nirx) sin (miry/b) . (4.115) 

m=l n=l 


The non-dimensional governing differential equation for the static problem takes 
the form 


00 oo 


2 ^ ^ = E L q sin (nirx) sin (miry/b) 

3x 9y rn=l n=l 


(4.116) 


where 


P 

%n - ■ - :2 - - • ^1 ^2 > 


4.5-1 Fully Clamped Plate 


The problem definition is completed by specifying the pertinent boundary con- 
ditions. As in the previous cases, Wq satisfies the zero deflection and zero moment 

conditions on all edges. The dynamic problem is the same as that presented in Sec- 
tion 4.3 and will not be repeated here. 

As outlined previously in Section 3, we write 


oo 



(4.117) 


for the outer solution as was done by Adeniji-Fashola and Oyediran [1] when they 
considered the same problem but for the fully-clamped rectangular plate case. The 
procedure followed is to solve equation (4.117) for v = 0,1,2,... subject to appro- 
priate boundary conditions. The fully-clamped case, which is discussed in Reference 
1 win not be repeated here. 
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4.5-2 Fully Simply- Supported Rectangvilar Plate 


For a fully simply- supported plate, it is easily seen that the boundary value 
problem for takes the form 


3 


2 

1 





sin (nirx) sin (muy/b) 


(4.118) 


Subject to i|)q vanishing on all boundaries. 

The solution takes the form 

sin (mrx) sin (miry/b) (4.119) 

where 


^mn = ^2'“'^'/^') • 

With relative ease, it can be verified that ip. must vanish while the solution to 
the t /)2 problem is 


ij '2 = fig sin (mrx) sin (miry/b) 


with 


«5 = 


'mn 


2 . ^2 2 /, 2.2 

(n 7T + m TT /b ) 


(3 


2 2 2 2 2 2 ,,. 

, n m^ n/h 


(4.120) 


This solution for the simply- supported case agrees with the solution by 
Timoshenko and Woinowski-Krieger [3]. 


4.5-3 Partly Clampled and Partly Hinged Plate [C-C-S-S] 

The governing differential equations are those presented in Reference 1 plus 
the pertinent boundary conditions. It is therefore clear that the solution takes the 
form 


i 

I ’^0 ~ ^mn (nirx) sin (miiy/b) 


(4.121) 
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’i'l ~ Lj'sinhaj^ sin (nirx) + L2' sinh (ct2x) sin (rmry/b) 


with 


(-1)™ Yn,n (imr/b) 

L^' = 521} 


V = 


^2 sinh 


(-1) riTT 

sinh «2 


(4 


(4 


The second order problem, ip2» can be shown to take the form 


.1. ,1, 

„ 2 “ '*'2 “ ’"2 4 


(4, 


subject to 


’p2(®»y) - " sin (oj^y) + L2' «2 (™ 7 ry/b) 

1 ri"^’ X 

’l'2(l»y) = “ [(“ 1 ) nTT Lj^' sinh(aj^y) - L2' 02 cosh ct2 sin (miry/b)] 

tl>2(x,0) - - [Lj' ctj^ sin (mrx) + (mir/b) L2' sinh (ct2x) 

2 

= - A m+1 

i|»2(x,l) ^2 [(mir/b) (- 1 ) L2’ sinh (a2x) - cosh (oj^b) sin (mrx) 


Clearly, we have the solution to as (see Appendix) 


il>2 - [Kj^ sinh (Hj^y) + K2 cosh (oTj^y) + sinh (aj^y) + E2 cosh (a^^y) 

+ Eg y cosh (otj^y)] sin (mrx) 

+ [Kg sinh (ogX) + cosh (otgX) + E^ sinh (otgX) + Eg cosh (ogX) 


. 122 ) 


.123) 


124) 


125) 
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+ Eg X cosh (a2x)] sin (imry/b) 


+ Kg sin (nirx) sin (miry/b) 


(4.126) 


The constants are determined in the Appendix and are listed in equation (A. 13). 


V. CONCLUDING REMARKS 


In this report, analytical results that include the small initial curvature effect 
on the free vibration of rectangular plates are presented in the limit when in- plane 
forces are much larger than the bending rigidity. These results are presented for 
various boundary conditions such as fully clamped, fully simply- supported, and 
partly clamped partly simply- supported cases. Different mathematical forms of the 
initial curvature are also considered. 

Taking advantage of the linear nature of the problem, a decomposition method 
is used to break the problem into two parts, viz: the dynamic and the static 

problems. The solution to the dynamic problem which is readily available in the 
literature is recovered independently by the present authors while the solutions to 
the static problem for various boundary conditions and initial curvature forms are 
presented. The method of matched asumptotic expansions is used. 

It is found that a small initial curvature has no effect on the natural frequency 
of these plates. However, the eigenmodes are greatly affected when compared with 
initially perfectly flat plates. The results presented here complement the earlier work 
of the present authors [ 1] . 
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APPENDIX 


We consider the fourth-order problem 


2 ^ ^^2 2 ^ ^^2 4 

i ^*^2 ^ ’ *0 


♦2<«>y) = ■ [^1 s™h (“jy) 


+ Lo a„ sin 


4 ^ 2 (l,y) = - g- (-l)’^ ^ niT sinh (a^y) - L2 02 cosh 


~ . / mTry \ 

“2 \ b / J 


ij/2(x,0) = 


ij;2(x,b) = 


^ |l. a-, sin (nir 

^2 L 


x) + ^ L2 sinh ( 


S2x)j 


(- 1 )*^^^ L2 sinh (a2x) - ctj^ cosh (oTj^b) sin ( 


\ 

mrx)| , 


where and L2 are known quantities (see text). 


To solve this boundary value problem, it is convenient to introduce the decom- 
position 

where satisfies the homogeneous part of equation (A.l) and the following boun- 

dary conditions: 


( 1 ) ~ 

'('2'' ''( 0 »y) = — — sinh (aj^y) 


sinh (otj^y) 
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sinh (agX) 

2 


1(^2^ ^\x,b) = ^ " g^b*^" ^2 («2^) 


(A. 4) 
(Concluded) 


( 2) 

The problem for 4>2 is described by 




3 ^, 1 . ( 2 ) 

2 ^ ’''2 _ „4 . 

2 - 2 ~ ^ '^O 

9y 


(A. 5) 


and the boundary conditions of equations (A. 2) less the portions satisfied by the 

f 2) 

problem in equations (A. 4). Thus, we have for the boundary conditions of the ip2 
problem 


( 1 ) 




^2 “ 


r 




^2 “2 


cosh a2 sin 


) 


<J»2^^\x,0) i sin (mrx) 




(A. 6) 


f ^1 ®1 

11^2'' ^'^(x,b) = — ^ cosh (aj^b) sin (mrx) 

2 


The boundary conditions dictate the solution of i/»2 


( 2 ) 


to take the form 




2<«(x.y) = 


n 


[Kj^ sinh (oTj^y) + K2 cosh (oj^y)] sin (mrx) + 2-# [Kg sin 


m 


(ot2x) 


+ cosh (o^gx)] 



+ K. 


sin (mrx) sin 



(A. 7) 
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with 


2L.| a. 

K, = — s — - coth (cTj^b) 




^3 = 


^2 


^1 

“1 

8 

2 

2 L2 

“2 

8 

1 

- 

“2 


coth a<j 


K. = 


/ 2 2 TT^\ 

Ymn ^ ^2 ) 


Bi 2 „2 ,2 * 


(A. 8) 


where y„ is defined in equation (4.119). 
mn ^ 

We attempt to solve for using the finite Fourier transform method. 

The transform with respect to x gives 


82^ ?2 ■ ^1^ h = ^1^ 'l^2^1.y) - <l^2^0,y)] , (A. 9) 

yy 

where iF2(n,y) is the transform of (|/2 with respect to x. Using equation (A. 4) in 
equation (A. 9), we have 


^2 ■ “1^ ^ 2 " ^ ^1 ^1 ^“1^^ • (A. 10) 

yy 

Thus, the solution to 'i'2^^^ form 

<l'2 (n,y) = Ej^ sinh (ctj^y) + E2 cosh (cTj^y) + Eg y cosh (a^y) (A. 11) 
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where compatibility requires that 


TT^ 

^2^ “l 


On taking the transform with respect to y, adding the result to equation (A. 11 ) 
and inverting, we obtain 


( 1 ) 


^2 (x,y) = sinh (a^y) + E2 cosh (a^^y) + y cosh (a^y)] sin (mrx) 


+ [E^ sinh (a2x) + Eg cosh (a2x) + Eg x cosh (a2x)] sin ( j (A. 12 ) 


( niTry ] 

V'^/ 


where 


E„ = 


2 2 - , 
rn 7 T p 2 -Li 2 

“2 


The unknowns Ej^, E2, E^, and Eg are, by requiring that ij)2^^\x,y) satisfy 
the boundary conditions of equation (A. 4 ), 

(-1)”“ mir L„ K ' L. a. 

E^ = = — — + ^ - — coth (S^b) - E„b coth (ar.b) 

^ Sgb sinh (a^b) ^^2 i 6 i 


®2 ~ ° 


(A. 13) 


(- 1 )” nir L2 ^2 


^4 = 


sinh 012 


23 i 


coth a„ - 


Eg coth «2 


where 


^2’ = 


(- 1 )''^^^ nir sinh 

~ 2 ^ „2 2 

a2 + n TT 


^4’ = 


(-l)m+l sinh 
~ 2 . m^ ir^ 

“1 ■'-72- 
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